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Abstract. A classical (or quantum) superintegrable system on an n-dimcnsional Rieman- 
nian manifold is an integrable Hamiltonian system with potential that admits 2n — 1 func- 
tionally independent constants of the motion that are polynomial in the momenta, the 
maximum number possible. If these constants of the motion arc all quadratic, the system is 
second order superintegrable. Such systems have remarkable properties. Typical properties 
are that 1) they are integrable in multiple ways and comparison of ways of integration leads 
to new facts about the systems, 2) they arc multiscparable, 3) the second order symmetries 
generate a closed quadratic algebra and in the quantum case the representation theory of the 
quadratic algebra yields important facts about the spectral resolution of the Schrodinger op- 
erator and the other symmetry operators, and 4) there are deep connections with expansion 
formulas relating classes of special functions and with the theory of Exact and Quasi-exactly 
Solvable systems. For n = 2 the author, E.G. Kalnins and J. Kress, have worked out the 
structure of these systems and classified all of the possible spaces and potentials. Here I 
discuss our recent work and announce new results for the much more difficult case n = 3. 
We consider classical superintegrable systems with nondegencratc potentials in three dimen- 
sions and on a conformally flat real or complex space. We show that there exists a standard 
structure for such systems, based on the algebra of 3 x 3 symmetric matrices, and that the 
quadratic algebra always closes at order 6. We describe the Stackel transformation, an in- 
vertible conformal mapping between superintegrable structures on distinct spaces, and give 
evidence indicating that all our superintegrable systems are Stackel transforms of systems 
on complex Euclidean space or the complex 3-sphere. We also indicate how to extend the 
classical 2D and 3D superintegrability theory to include the operator (quantum) case. 
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1 Introduction and examples 

In this paper I will report on recent and ongoing work with E.G. Kalnins and J. Kress to uncover 
the structure of second order superintegrable systems, both classical and quantum mechanical. 
I will concentrate on the basic ideas; the details of the proofs can be found elsewhere. The 
results on the quadratic algebra structure of 3D conformally flat systems with nondegenerate 
potential have appeared recently. The results on the 3D Stackel transform and multiseparability 
of superintegrable systems with nondegenerate potentials are announced here. 

Superintegrable systems can lay claim to be the most symmetric solvable systems in mathe- 
matics though, technically, many such systems admit no group symmetry. In this paper I will 
only consider superintegrable systems on complex conformally flat spaces. This is no restriction 
at all in two dimensions. An n-dimensional complex Riemannian space is conformally flat if 
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and only if it admits a set of local coordinates xi,. . . ,Xn such that the contravariant metric 

n 

tensor takes the form g^^ = 5^^ /\{x). Thus the metric is ds^ = A(a;)( ^ dxj). A classical 

1=1 

superintegrable system H = ^ g^^'PiVj + ^{'^) o^i the phase space of this manifold is one that 

admits 2n — 1 functionally independent generalized symmetries (or constants of the motion) 5^, 
k = 1, . . . , 2n — 1 with Si = 7i where the Sk are polynomials in the momenta pj. That is, 
{H, Sk} = where 

n 

is the Poisson bracket for functions f{x,p), g{x,p) on phase space ^|HI31l3ISliniIZllHl- It is 

easy to see that 2n — l is the maximum possible number of functionally independent symmetries 
and, locally, such (in general nonpolynomial) symmetries always exist. The system is second 
order superintegrable if the 2n — 1 functionally independent symmetries can be chosen to be 
quadratic in the momenta. Usually a superintegrable system is also required to be integrable, 
i.e., it is assumed that n of the constants of the motion are in involution, although I will not 
make that assumption in this paper. Sophisticated tools such as i?-matrix theory can be applied 
to the general study of superintegrable systems, e.g., [HI El ^5- However, the most detailed 
and complete results are known for second order superintegrable systems because separation 
of variables methods for the associated Hamilton-Jacobi equations can be applied. Standard 
orthogonal separation of variables techniques are associated with second-order symmetries, e.g., 
jl2l 1131 1141 1151 1161 117j and multiseparable Hamiltonian systems provide numerous examples of 
superintegrability. Thus here I concentrate on second-order superintegrable systems, on those 
in which the symmetries take the form S = '^o}^ {x)piPj + W{x), quadratic in the momenta. 

There is an analogous definition for second-order quantum superintegrable systems with 
Schrodinger operator 

H = ^ + V{x), A = ^ J] 5., (V59*^)5,, , 

the Laplace-Beltrami operator plus a potential function J^l- Here there are 2n — 1 second-order 
symmetry operators 

Sk = ^Yl + W^'^Hx), fc = 1, . . . , 2n - 1 

with Si = H and [H, S^] = HS^ — S^H = 0. Again multiseparable systems yield many 
examples of superintegrability, though not all multiseparable systems are superintegrable and 
not all second-order superintegrable systems are multiseparable. 

The basic motivation for studying superintegrable systems is that they can be solved explicitly 
and in multiple ways. It is the information gleaned from comparing the distinct solutions and 
expressing one solution set in terms of another that is a primary reason for their interest. 

Two dimensional second order superintegrable systems have been studied and classified by 
the author and his collaborators in a recent series of papers [1811191 QUI 121j . Here we concentrate 
on three dimensional (3D) systems where new complications arise. We start with some simple 
3D examples to illustrate some of the main features of superintegrable systems. (To make clearer 
the connection with quantum theory and Hilbert space methods we shall, for these examples 
alone, adopt standard physical normalizations, such as using the factor — ^ in front of the free 
Hamiltonian.) Consider the Schrodinger equation = oi {h = m = 1, xi = x, X2 = y, 
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The generahzed anisotropic oscillator corresponds to the 4-parameter potential 



^2 

Vix, y,z) = — {x^ + y2 + 4(z + pf) + ^ 
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(This potential is "nondegenerate" in a precise sense that I will explain later.) The corre- 
sponding Schrodinger equation has separable solutions in five coordinate systems: Cartesian 
coordinates, cylindrical polar coordinates, cylindrical elliptic coordinates, cylindrical parabolic 
coordinates and parabolic coordinates. The energy eigenstates for this equation are degenerate 
and important special function identities arise by expanding one basis of separable eigenfunc- 
tions in terms of another. A second order symmetry operator for this equation is a second order 
linear differential operator S such that [H, S] = 0, where [A, B] = AB — BA. A basis for these 
operators is 



i-2 _ 1 

71^ o2 2 2,1 4 

Mi = d^- UJ X H 



Mo 



a2 2 2 

dy-oj y 



^^2 4 



P = dl-^uj'^{z + pf, 



T - T"^ 



^2 



1\ x' 



Si 



2 2 



-idxLi'i + Li39^.) + pdi + (z + p) uj'x 



ft-l 4 



^2 

where L. 



--{dyL23 + L2zdy) + pdl + {z + p)[ u^y^ 



'^idxi — Xjdx^- It can be verified that these symmetries generate a "quadratic 



algebra" that closes at level six. Indeed, the nonzero commutators of the above basis are 

[Mi,L] = [L,M2] =Q, [L,Si] = [S2,L]=B, [M,,S^] = Ai, [P,S^] = -A,. 

Nonzero commutators of the basis symmetries with Q (4th order symmetries) are expressible in 
terms of the second order symmetries: 

[M„ Q] = [Q, M2] = 4{Mi, M2} + IQuj^L, [Si, Q] = [Q, S^] = 4{Mi, M2}, 
[L, Q] = 4{Mi, L} - 4{M2, L] + 16 (l - kl) Mi - 16 (l - kl) M2. 

There are similar expressions for commutators with B and the Ai. Also the squares of Q, B, 
Ai and products such as {Q,B}, (all 6th order symmetries) are all expressible in terms of 2nd 
order symmetries. Indeed 



Q2 = _{L, Ml, M2} + 8uj'^{L, L} - 16 (1 - kf) - 16 (l - k'^) M| 
3 

+ y{Mi, M2} - ^iv'L - 128^2 (1 _ ^2) (1 _ ^2) ^ 

{Q, B} = -^{M2, L, Si} - ^{Mi, L, ^2} + 16 (1 - kj) {M2, 52} + 16 (1 - k^) {Mi, Si} 

64 64 
--{Mi,52}--{M2,5i}. 
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Here {Ci, ... ,Cj} is the completely symmetrized product of operators Ci, . . . ,Cj. (For complete 
details see [H].) The point is that the algebra generated by products and commutators of the 
2nd order symmetries closes at order 6. This is a remarkable fact, and ordinarily not the case 
for an integrable system. 

A counterexample to the existence of a quadratic algebra in Euclidean space is given by the 
Schrodinger equation with 3-parameter extended Kepler-Coulomb potential: 



dx"^ dy'^ dz^ 



2a 

2E + 




^' = 0. 



This equation admits separable solutions in the four coordinates systems: spherical, sphero- 
conical, prolate spheroidal and parabolic coordinates. Again the bound states are degenerate 
and important special function identities arise by expanding one basis of separable eigenfunctions 
in terms of another. However, the space of second order symmetries is only 5 dimensional and, 
although there are useful identities among the generators and commutators that enable one to 
derive spectral properties algebraically, there is no finite quadratic algebra structure. The key 
difference with our first example is, as we shall show later, that the 3-parameter Kepler-Coulomb 
potential is degenerate and it cannot be extended to a 4-parameter potential. 

In |2fl[ l21j there are examples of superintegrable systems on the 3-sphere that admit a quadra- 
tic algebra structure. A more general set of examples arises from a space with metric 

ds^ = X{A, B, C, D, x) {dx^ + dy'^ + dz^) , 

where 

/3 z\ f 1 3z 

\ = A{x + iy) + B{-{x + iyf + -\ + C {{x + iyf + — (x - iy) + — (x + iy) 



The nondegenerate classical potential \s V = A(a, /3, 7, 6, x)/\{A, B,C,D,x). A = B = C = 
D = this is a nondegenerate metric on complex Euclidean space. The quadratic algebra always 
closes, and for general values of A, B, C, D the space is not of constant curvature. As will be 
apparent later. This is an example of a superintegrable system that is Stackel equivalent to 
a system on complex Euclidean space. 

Observed common features of superintegrable systems are that they are usually multisepara- 
ble and that the eigenfunctions of one separable system can be expanded in terms of the eigen- 
functions of another. This is the source of nontrivial special function expansion theorems |23j . 
The symmetry operators are in formal self-adjoint form and suitable for spectral analysis. Also, 
the quadratic algebra identities allow us to relate eigenbases and eigenvalues of one symmetry 
operator to those of another. The representation theory of the abstract quadratic algebra can be 
used to derive spectral properties of the second order generators in a manner analogous to the 
use of Lie algebra representation theory to derive spectral properties of quantum systems that 
admit Lie symmetry algebras, |23 | I24[ l25l I26j . (Note however that for superintegrable systems 
with nondegenerate potential, there is no first order Lie symmetry.) 

Another common feature of quantum superintegrable systems is that they can be modified by 
a gauge transformation so that the Schrodinger and symmetry operators are acting on a space of 
polynomials |22j. This is closely related to the theory of exactly and quasi-exactly solvable sys- 
tems |281 129| . The characterization of ODE quasi-exactly solvable systems as embedded in PDE 
superintegrable systems provides considerable insight into the nature of these phenomena j3Uj . 

The classical analogs of the above examples are obtained by the replacements dx^ — > Pxi and 
modification of the potential by curvature terms. Commutators go over to Poisson brackets. 
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The operator symmetries become second order constants of the motion. Symmetrized operators 
become products of functions. The quadratic algebra relations simplify: the highest order terms 
agree with the operator case but there are fewer nonzero lower order terms. 

Many examples of 3D superintegrable systems are known, although they have not been clas- 
sified 1^ EH ESI EH EHl EH] ■ Here, we employ a theoretical method based on integrability 
conditions to derive structure common to all such systems, with a view to complete classifica- 
tion, at least for classical systems with nondegenerate potentials. We show that for systems 
with nondegenerate potentials there exists a standard structure based on the algebra of 3 x 3 
symmetric matrices, and that the quadratic algebra closes at level 6. For 2D nondegenerate 
superintegrable systems we earlier showed that the 3 = 2(2) — 1 functionally independent con- 
stants of the motion were (with one exception) also linearly independent, so at each regular 
point we could find a unique constant of the motion that matches a quadratic expression in 
the momenta at that point. However, for 3D systems we have only 5 = 2(3) — 1 functionally 
independent constants of the motion and the quadratic forms span a 6 dimensional space. This 
is a major problem. However, for nondegenerate potentials we prove the "5 implies 6 Theorem" 
to show that the space of second order constants of the motion is in fact 6 dimensional: there is 
a symmetry that is functionally dependent on the symmetries that arise from superintegrability, 
but linearly independent of them. With that result established, the treatment of the 3D case 
can proceed in analogy with the nondegenerate 2D case treated in ^Sj. Though the details are 
quite complicated, the spaces of truly 2nd, 3rd, 4th and 6th order constants of the motion can 
be shown to be of dimension 6, 4, 21 and 56, respectively and we can construct explicit bases 
for the 4th and 6th order constants in terms of products of the 2nd order constants. This means 
that there is a quadratic algebra structure. 

Using this structure we can show that all 3D superintegrable systems with nondegenerate 
potential are multiseparable. We study the Stackel transform, or coupling constant metamor- 
phosis [371 138j . for 3D classical superintegrable systems. This is a conformal transformation of 
a superintegrable system on one space to a superintegrable system on another space. We give 
evidence that all nondegenerate 3D superintegrable systems are Stackel transforms of constant 
curvature systems, just as in the 2D case, though we don't completely settle the issue. This 
provides the theoretical basis for a complete classification of 3D superintegrable systems with 
nondegenerate potential, a program that is underway. Finally we indicate the quantum analogs 
of our results for 3D classical systems. 

2 Conformally flat spaces in three dimensions 

We assume that there is a coordinate system x, y, z and a nonzero function \{x^ y, z) = 
exp G{x, y, z) such that the Hamiltonian is 




A quadratic constant of the motion (or generalized symmetry) 



3 



S=Y^ a^Hx, y, z)pkPj + W{x, y,z) = C + W, 



a- 



= a 



kj 



k,j=l 



must satisfy {7i,S} = 0, i.e. 



^ 7^ J, 



6 



W. Miller 



+ Oj* + = 0, i, j, k distinct 

and 

3 

Wk = \Y,a'^ys. A; = 1,2, 3. (1) 

s=l 

(Here a subscript j denotes differentiation with respect to Xj.) The requirement that dxfWj = 
dxjWi, i ^ j leads from to the second order Bertrand-Darboux partial differential equations 
for the potential. 



s=l 



VsjXa'' - VsAa'' + Vs (^{Xa'')j - (Aa^^),)] = 0. (2) 



For second order superintegrabilty in 3D there must be five functionally independent con- 
stants of the motion (including the Hamiltonian itself). Thus the Hamilton-Jacobi equation 
admits four additional constants of the motion: 

3 

5h = ^ alhfkPj + W^h) = + W(^h), h = l,...,4. 
j,k=i 

We assume that the four functions together with TC are functionally independent in the six- 
dimensional phase space. (Here the possible V will always be assumed to form a vector space 
and we require functional independence for each such V and the associated W^^\ This means 
that we require that the five quadratic forms "Hq are functionally independent.) In [20] it 
is shown that the matrix of the 15 Bertrand-Darboux equations for the potential has rank at 
least 5, hence we can solve for the second derivatives of the potential in the form 

V22 = Vii + A^^Vi + B^Vs + C72V3, 
V33 = Vii + A^^Vi + B^^V2 + ^33^3, 
V12 = A^Vi + B^Vs + C^Va, 
^13 = A^^Vi + B^^V2 + C^^Vs, 

V23 = A^^Vi + B^^V2 + ^23^3. (3) 

If the matrix has rank > 5 then there will be additional conditions of the form Dj^^-^Vi + D'^^^V2 + 
D^s)^3 = 0. Here the A^^ , B^^ , C*-^, 1?^^^ are functions of x that can be calculated explicitly. For 
convenience we take A'^ = A^\ B'i = B^\ C'^ = 0\ 

Suppose now that the superintegrable system is such that the rank is exactly 5 so that the 
relations are only (jSJ. Further, suppose the integrability conditions for system (jSJ are satisfied 
identically. In this case we say that the potential is nondegenerate. Otherwise the potential is 
degenerate. If V is nondegenerate then at any point xq, where the A'^^ , B^^ , C*-' are defined and 
analytic, there is a unique solution V{x) with arbitrarily prescribed values of Vi{xq), V2{xq), 
V3{xq), Vii{xq) (as well as the value of V{xo) itself.) The points a^o are called regular. The points 
of singularity for the A^^ , B^^ , C^^ form a manifold of dimension < 3. Degenerate potentials 
depend on fewer parameters. For example, it may be that the rank of the Bertrand-Darboux 
equations is exactly 5 but the integrability conditions are not satisfied identically. This occurs 
for the generalized Kepler-Coulomb potential. 

From this point on we assume that V is nondegenerate. Substituting the requirement for 
a nondegenerate potential © into the Bertrand-Darboux equations ((2)) we obtain three equa- 
tions for the derivatives , the first of which is 



- ) V, + (42 - af ) V2 + (43 - a33)V3 
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- (A^Vi + S^Va + C^Vg) + (^33^^ ^ ^33y^ ^ ^33y^^ 

= (-G30" + Gioi3) + (-Gga^^ + Gia23) ^3 + (-Gsa^^ + da^^) V3, 

and the other two are obtained in a similar fashion. 

Since V is a nondegenerate potential we can equate coefficients of Vi, V2, V3, Vn on each side 
of the conditions 51^23 = c^aVis = 83X^2, 193V23 = 82^33, etc., to obtain integrability conditions, 
the simplest of which include 

A^-^=B^^ = C^^ 5l2_^22^(^13_^33^ 
523 ^ ^31 ^ ^22^ ^23 ^ ^12 ^ ^33^ 

+ B^^A^^ + Af + A^^A^^ + B^^A^^ + C^^A^^ = Af + B^^A^^ + C^^A^^, 
Al^ + A^^A^^ + B^'-^A^^ + C^^A^^ = Af + ^23^12 ^ (^23^13 

= ^12 ^ ^13^12 ^ ^12^23 ^ ^12^33_ 

Using the nondegenerate potential condition and the Bertrand-Darboux equations we can 
solve for all of the first partial derivatives aj^ of a quadratic symmetry to obtain 

(4) 





-Gia^i - 


G2a^2 _ ^^^13^ 


22 
a2 = 


-Giai2 - 


G2a22 - G3a2^ 


„33 _ 
"3 — 


-Giai3 - 


G2023 - Gaa^s, 


3af = 


= ai2A22 - 


(a22-aii)Ai2-a23Ai3^ai3^23 




+ Gaa" 





3a^i = -2a'^A^^ + 2 (a^^ - a") + 2a23Ai3 _ 2^13^23 
-2G2aii + Giai2-G2a22-G3a23, 

= -0^2^23 ^ (^33 _ ^11) ^13 ^ ^23^12 _ ^13^33 

- Gia" - Gaa^^ - 2030^3 ^ ,^^^33^ 

3af = 2a^2(^23 _ 2 (^33 _ ^11^ ^13 _ 2a^'^C^^ + 2a^^C^^ 

- Gia" - G2a^^ + Gsa^^ - 2Gia^^, 

3af = a23(i?33 _ ^22) _ (^33 _ ^22) ^23 _ ^13^12 ^ ^12^13 

- Gia^^ - 2G2a^^ - ^30^^ + G3a^^ 

3af = -20^^ {B^^ - 5^2) + 2(a^2 - a^^) + 20^25^2 _ 2ai2si3 

- Gia^^ + ^20^3 - Gsa^^ - 2G3a^^ 

3af = -a23Ai2 + (a^i - a33) A^'^ + ai3A33 + ai2^23 

- 2Gia^2 - G2a23 - Gsa^^ + G3a^\ 

3aii = 2a23yli2 + 2 (a^s - a'') A^^ - 2a'' A^^ - 2a'^A''' 

+ Gia^' - G2a^' - Gsa'^ - 2G3a'\ 
3af = -2ai3ci2 + 2 (a^^ - a'^) + 2ai2C7i3 - 2a23 (C722 _ ^33) 

- Gia'^ - G2a^^ + Gsa'^' - 2G2a^^ 

3af = ai3ci2 - (a22 - a33) G^' - a'^G^' - (G^^ - G^^) 

- Gia'^ - G2a^^ - 2^30^^ + G2a^^ 

342 = -al3523 + („22 _ ^11) ^12 _ ^12^22 ^ ^23^13 
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- da" - 2G2a^^ - Gga^^ ^ q^^22 ^ 

3af = 2a^^B'^^ - 2 (a^^ - a") _^ 2a^2s22 - 2a2^5^^ 

- Gia" + G2a^2 _ (^^^^is _ 2^10^2, 

3af = al2 (i?23 + ^22) ^ ^11 (^13 ^ ^12) _ ^22^12 _ ^33^13 

+ a^3 (^33 + C723) - (C713 ^ ^12^ _ 2Gia23 + Gza^^ + Gga^^^ 

3aJ2 = ^12 (_2ij23 ^ ^22) ^ ^11 (^12 _ 2^13) _ ^22^12 ^ 2a33i?13 

+ a^^ (-2^33 + C723) + (-C^3 ^ 2^12^ _ 2030^^ + €20^^ + da^^. 

3a^3 = ^12 (^23 _ 2^22) ^ ^11 (^13 _ 2^12) + 2022^^2 _ ^33^13 

+ a^3 (^33 _ 2(^23^ ^ ^23 (^2(^13 _ ^12^ _ 2G2a^3 ^ (^^^23 ^ ,^^^12^ 

plus the linear relations 

^23 ^ ^13 ^ (^12^ ^23 _ ^31 _ ,^22 ^ 

^12 _ ^22 ^ ^33 _ (jVi ^ ^33 ^ ^12 _ ^23 ^ g 

Using the linear relations we can express C^"^, C^^, C'^'^, C'^^ and B^^ in terms of the remaining 
10 functions. 

Since the above system of first order partial differential equations is involutive the general 
solution for the 6 functions a^'' can depend on at most 6 parameters, the values a^^{xQ) at a fixed 
regular point Xq. For the integrability conditions we define the vector- valued function 

h{x,y,z) = (all ^12 ^13 ^22 ^23 ^33) 

and directly compute the 6x6 matrix functions A^^^ to get the first-order system 

d^.h = A^^^h, j = 1,2,3. 
The integrability conditions for this system are are 

J^^h - Afh = A^'^A'^^^h - A'^^^A'^^h = (5) 
In terms of the 6x6 matrices 

5(1) = - - [A^'\A^% 5(2) = 41) - _ [^(3)^_4(i)]^ 

5(3) =42)- 41) -[^(1), ^(2)], 

the integrabilty conditions are 

5(i)/i = 5(2)/i = 5(3)/i = 0. (6) 



3 The 5 =^ 6 Theorem 

Now assume that the system of equations @ admits a 6-parameter family of solutions a^''. (The 
requirement of superintegrability appears to guarantee only a 5-parameter family of solutions.) 
Thus at any regular point we can prescribe the values of the a^'' arbitrarily. This means that © 
or Q holds identically in h. Thus 5(1) = 5(2) = 5(3) = 0. This would be the analog of 
what happens in the 2D case where there are Sindependent terms in the quadratic form and 
3 functionally (and linearly) independent symmetries. However, in the 3D case there are only 
5 functionally independent symmetries, so we can't guarantee that the symmetry equations 
admit a 6-parameter family of solutions. Fortunately, by careful study of the integrability 
conditions of these equations and use of the requirement that the potential is nondegenerate, 
we can prove the 5 6 theorem PU] . 
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Theorem 1 (5 =^ 6). Let V be a nondegenerate potential corresponding to a conformally 
flat space in 3 dimensions that is superintegrable, i.e., suppose V satisfies the equations Q 
whose integrability conditions hold identically, and there are 5 functionally independent con- 
stants of the motion. Then the space of second order symmetries for the Hamiltonian 7i = 
{Px~'rPy+p1) /^{x,y,z) + V{x,y,z) (excluding multiplication by a constant) is of dimension 
D = 6. 

Corollary 1. IfTi + V is a superintegrable conformally flat system with nondegenerate potential, 
then the dimension of the space of 2nd order symmetries 

3 

5 = ^ a''\x,y,z)pkPj + W{x,y,z) 
k,j=i 

is 6. At any regular point {xo,yQ, zq), and given constants a^^ = a^^ , there is exactly one 
symmetry S (up to an additive constant) such that a^^ {xo,yQ, zq) = a^^ . Given a set of 5 
functionally independent 2nd order symmetries C = {Sg : i = 1, ... 5} associated with the 
potential, there is always a 6th second order symmetry Sq that is functionally dependent on C, 
but linearly independent. 

4 Third order constants of the motion 

The key to understanding the structure of the space of constants of the motion for superintegrable 
systems with nondegenerate potential is an investigation of third order constants of the motion. 
We have 

3 

/C = ^ a''^\x, y, z)pkPjPi + h^{x, y, z)pi, 

k,j,i=l 

which must satisfy {Ti,JC} = 0. Here a^^^ is symmetric in the indices k, j, i. 
The conditions are 

s 

3af + af = -3 ^a^'^' (In A)„ i^j 

s 

+ af = -i 5] a^^^ (In A). - ^ a-^(ln A)„ i^ j, 

s s 

2af *^ + af + af = ^)^' ^' ^ distinct, 

s 

bi + b'j=3XYa"'^Vs, j^k, j,fc = 1,2,3, 

s 

6^: = ^ A E a'^^V, - i E 6^ (In A)„ j = 1, 2, 3, 

s s 

and 

Yb'Vs = 0. 

s 

The a^^^ is just a third order Killing tensor. We are interested in such third order symmetries 
that could possibly arise as commutators of second order symmetries. Thus we require that 
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the highest order terms, the af'^^ in the constant of the motion, be independent of the four 
independent parameters in V . However, the must depend on these parameters. We set 

3 

b\x,y,z) = ^f'^{x,y,z)Vj{x,y,z). 
i=i 

(Here we are excluding the purely first order symmetries.) In j2Uj the following result is obtained. 

Theorem 2. Let K, he a third order constant of the motion for a conformally flat superintegrable 
system with nondegenerate potential V : 

3 3 
/C = ^ a''^'\x, y, z)pkPjPi + ^ b\x, y, z)pi. 

k,j,i=l 1=1 

Then 

3 

b\x,y,z) = ^f'^'^{x,y,z)Vj{x,y,z) 
i=i 

with Z^'-' + /■''^ = 0, 1 < £, j < 3. The a^^^ , are uniquely determined by the four numbers 

f^''^{xo,yo,zo), f^'^{xo,yo,zo), /^'^(xq, yo, ^o), /^'^(a^o, yo, ^^o) 
at any regular point (xq, yo,zo) ofV. 
Let 

cSi = ^ a'^l-^PkPj + W(i) , '^2 = ^ <\i)PkPj + W{2) 

be second order constants of the the motion for a superintegrable system with nondegenerate 
potential and let A{i){x,y, z) = {_a^^i^{x^y^ z)^ , i = 1,2 be 3 x 3 matrix functions. Then the 
Poisson bracket of these symmetries is given by 

3 

{81,82}= ^ a^^\x,y,z)pkPjPi + b\x,y,z)pe, 

k,j,i=l 

where 

fk,l _ X (kj jl kj jl \ 

J - 2A 2^ l«(2)"(l) - "(l)"(2) j • 

3 

Differentiating, we find 

/f = 2A ^ (^iaf2)«(i) + ~ ^i"a)«(2) ~ "a)^^"(2)) + G^f''''■. (7) 

3 

Clearly, {5i,52} is uniquely determined by the skew-symmetric matrix 

[^(2),^(l)] = ^(2)^(1) - ^(1)^(2), 

hence by the constant matrix [A{^2){xQ-,yQ-, zq) , A(^i){xQ,yQ, zq)] evaluated at a regular point, and 
by the number T{xq, yo, Z)) = fl'^{xo,yo, zq). 

For superintegrable nondegenerate potentials there is a standard structure allowing the iden- 
tification of the space of second order constants of the motion with the space of 3 x 3 
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symmetric matrices, as well as identification of the space of third order constants of the mo- 
tion with a subspace of the space x F of 3 x 3 skew-symmetric matrices crossed with 
the line F = {J^(xo)}. Indeed, if xq is a regular point then there is a 1 — 1 linear correspon- 
dence between second order symmetries S and their associated symmetric matrices ^(a^o)- Let 
{61,62}' = {62,61} be the reversed Poisson bracket. Then the map 

{6^,62}' ^[A^i){xo),A^2){xo)] 

is an algebraic homomorphism. Here, 61, 62 are in involution if and only if matrices A(^i^{xq), 
A[2){^o) commute and J-'{xq) = 0. If {5i,52} 7^ then it is a third order symmetry and can 
be uniquely associated with the skew-symmetric matrix [A(^i-^{xQ),A(^2)i^o)] the parame- 
ter J-{x()) . Let be the 3x3 matrix with a 1 in row i, column j and for every other matrix 
element. Then the matrices 

j^iij) ^ ^^gij + gji-^ ^ ^00^ ij = 1, 2, 3 (8) 
form a basis for the 6-dimensional space of symmetric matrices. Moreover, 

where 

j^iij) ^ _ gji^ = i,j = 1,2, 3. 

Here B^''^ = and B^^^\ B^^^\ ^(^i) fo^m a basis for the space of skew-symmetric matrices. To 
obtain the commutation relations for the second order symmetries we need to use relations ^ 
to compute the parameter J-'{xq) associated with each commutator [^(^-'^ ^('^^^j. The results 
are straightforward to compute, using relations @. 



Commutator 




3^/A 






,^(")] 


= fi(21) 


_3^13 _ ^23 _ ^3 




[^(13) 


,^(")] 


= 5(31) 


^12 _ 533 + 




[^(22) 


,^(")] 


= 


-4yl23 




[^(23) 


,^(")] 


= 


2(>122 _ ^33) 




[^(33) 


,^(")] 


= 


4^23 








= i5(32) 


i(3Sl2_^22^3^33_ 


Gi) 


[^(22) 


,^(12)] 


= 5(21) 


-3^23 _ ^13 _ 




[^(23) 


,^(12)] 


= i5(31) 


1(_3533_3^12^2522^(^2) 


[^(33) 


,^(12)] 


= 


2(^23 - yll3) 




[^(22) 


,^(13)] 


= 


-2533 




[^(23) 


,^(13)] 


= i5(21) 


_C33 + 1^23 _ 1^13 _ 




[^(33) 




= 5(31) 


^12 ^ ^33 ^ 




[^(23) 


,^(22)] 


= ^(32) 


^33 _ ^22 _ ^12 _ 




[^(33) 


,^(22)] 


= 


-4^23 




[^(33) 


,^(23)] 


= ^(32) 


^22 _ ^33 _ 512 _ 





A consequence of these results is (20] 
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Corollary 2. Let V he a superintegrahle nondegenerate potential on a conformally flat space, 
not a Stdckel transform of the isotropic oscillator. Then the space of truly third order constants 
of the motion is 4- dimensional and is spanned by Poisson brackets of the second order constants 
of the motion. 

Corollary 3. We can define a standard set of 6 second order basis symmetries 



5 Maximum dimensions of the spaces of polynomial constants 

In order to demonstrate the existence and structure of quadratic algebras for 3D nondegenerate 
superintegrable systems on conformally flat spaces, it is important to compute the dimensions 
of the spaces of symmetries of these systems that are of orders 4 and 6. These symmetries are 
necessarily of a special type. The highest order terms in the momenta are independent of the 
parameters in the potential, while the terms of order 2 less in the momenta are linear in these 
parameters, those of order 4 less are quadratic, and those of order 6 less are cubic. We will 
obtain these dimensions exactly, but first we need to establish sharp upper bounds. 

The following results are obtained by a careful study of the defining conditions and the 
integrability conditions for higher order constants of the motion [20] : 

Theorem 3. The maximum possible dimension of the space of purely fourth order symmetries 
for a nondegenerate 3D potential is 21. The maximal possible dimension of the space of truly 
sixth order symmetries is 56. 

6 Bases for the fourth and sixth order constants of the motion 

It follows from Section [S] that, for a superintegrable system with nondegenerate potential, the 
dimension of the space of truly fourth order constants of the motion is at most 21. Note from 
Section |l]that at any regular point Xo, we can define a standard basis of 6 second order constants 
of the motion S^^^^ = A^^^^ + W^^^^ where the quadratic form A^^^^ has matrix A^'^^^ defined by ^ 
and W^^^^ is the potential term with W^^^\xq) = identically in the parameters W^°'\ By taking 
homogeneous polynomials of order two in the standard basis symmetries we can construct fourth 
order symmetries. 

Theorem 4. The 21 distinct standard monomials S^^^^S^^''\ defined with respect to a regular 
point xq, form a basis for the space of fourth order symmetries. 

Indeed, we can choose the basis symmetries in the form 





for i, j,k = 1 



. . . ,3, i, j , k pairwise distinct (15 possibilities). 




for i,j 



3 



1, . . . ,3, i, j pairwise distinct (3 possibilities). 



for i,j, 



k = 1 



. . . ,3, i, j , k pairwise distinct (3 possibilities). 
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It is a straightforward computation to show that these 21 symmetries are hnearly independent. 
Since the maximum possible dimension of the space of fourth order symmetries is 21, they must 
form a basis. See j^HI for the details of the proof. 

Now from Section |21 the dimension of the space of purely sixth order constants of the motion 
is at most 56. Again we can show that the 56 independent homogeneous third order polynomials 
in the symmetries S^^^^ form a basis for this space. 

At the sixth order level we have the symmetries 

1. (5^"))'^, (cS^**))^^^-'-'), (^(^O^^^O'fc) 

for i,j,k = l,...,3, i, j, k pairwise distinct (18 possibilities). 

for i, J, A; = 1, . . . , 3, i, j, k pairwise distinct (10 possibilities). 

3. 5(^™)(cS(^^)5(^'^') - {S'-'^^f) 

for i, j = 1, . . . ,3, i, j pairwise distinct (10 possibilities). 

4. (5fe')5(*'=) _ ^(^O^O'fc)) 

for i,j,k = 1, . . . , 3, i, j, k pairwise distinct (18 possibilities). 

Theorem 5. The 56 distinct standard monomials S^^^^S^^^^S^^"^\ defined with respect to a re- 
gular xq, form a basis for the space of sixth order symmetries. 

See pUj for the details of the proof. We conclude that the quadratic algebra closes. 



7 Second order conformal Killing tensors 

There is a close relationship between the second-order Killing tensors of a conformally flat space 
in 3D and the second order conformal Killing tensors of flat space. A second order conformal 
Killing tensor for a space V with metric ds"^ = X{x)(^dx1 + dx^ + dx"^) and free Hamiltonian 
= {Pi + P2 +^1)/^ is a quadratic form S = J2 i^i^ ^2, X3)piPj such that 

{n,S} = f{xi,X2,X3)n, 

form some function /. Since / is arbitrary, it is easy to see that 5 is a conformal Killing tensor 
for V if and only if it is a conformal Killing tensor for flat space dxf + + The conformal 
Killing tensors for flat space are very well known, e.g., JH]- The space of conformal Killing 
tensors is infinite dimensional. It is spanned by products of the conformal Killing vectors 

Pi, P2, P3, X3P2-X2P3, X1P3-X3P1, X2P1-X1P2, X1P1+ X2P2 + X3P3, 
[xl - x\- x|) pi + 2X1X3^3 + 2X1X2P2, {x\ - x\- x|) P2 + 2X2^3^3 + 2x2XiPi, 
{xl - xl- X2) P3 + 2X3X1^1 + 2X3X2P2, 

and terms ^(xi, X2, X3) (pf +P2 +P3) where g is an arbitrary function. Since every Killing tensor 
is also a conformal Killing tensor, we see that every second-order Killing tensor for V3 can be 
expressed as a linear combination of these second-order generating elements though, of course, 
the space of Killing tensors is only finite dimensional. This shows in particular that every a^^ 
and every a** — a^^ with i 7^ j is a polynomial of order at most 4 in xi, X2, X3, no matter what 
is the choice of A. 

It is useful to pass to new variables 

a^' 

for the Killing tensor, where a'^^ = q}'^ — a^^, a^^ = a^^ — a^^. Then we can establish the result 
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Theorem 6. Necessary and sufficient conditions that the quadratic form 5 = ^'^''PiPj 

a second order Killing tensor for the space with metric ds"^ = \{^dx\ + dx"^ + c?x|) are 

1. S is a conformal Killing tensor on the flat space with metric dx\ + dx"^ + dx\. 

2. The following integrability conditions hold: 

(Aaa^^ + X3a^^)^ = {Xia^^ + {a'^^X)^ + ha^^)^ , 

(Aiai2 + (a^^A)^ + Asa^^)^ = {X,a'^ + X.a''^ + (a^^'X)^)^ . 

8 The Stackel transform for three-dimensional systems 

The Stackel transform j^Z] or coupling constant metamorphosis [SHj plays a fundamental role 
in relating superintegrable systems on different manifolds. Suppose we have a superintegrable 
system 

X{x,y,z) 

in local orthogonal coordinates, with nondegenerate potential V{x,y,z): 

V33 = Vii + A^^Vi + B^^V2 + C^^Vs, 
V22 = Vu + A^Vi + ^22^2 + C72V3, 
V23 = A^^Vi + ^23^2 + C'^^V^, 
Fi3 = A^^Vi + B^^V2 + C^Va, 

F12 = ^^Vi + S^Va + C^Va (9) 

and suppose U{x,y,z) is a particular solution of equations ©, nonzero in an open set. Then 
the transformed system 

X{x,y,z) 

with nondegenerate potential V{x,y,z): 

V33 = Vii + A^^Vi + B^^V2 + C^^Vs, 
V22 = Vii + A^'^Vi + ^22-(/2 + (522^3, 
1/23 = A^'^Vi + ^23^2 + C^^Vs, 

F12 = i^^Vi + B^^V2 + C^^Va, 
is also superintegrable, where 

V 

A = AC/, V = -, 

i33^^33^2^^ ^33^^33^ ^33 = ^33- 2^, 

i22^^22^2^^ ^22^^22_2^^ ^22 = ^22, 
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^23 


= A^\ 


^23 


= ^23- 


U3 

U ' 


^23 


= C23 


^13 


= A^^- 


t/3 
[/' 






^13 


= C^3 


^12 


= A^^- 


U2 
[/' 











U2 

" u' 

Let S* = a^^piPj + W = Sq + W he a second order symmetry of H and 81/ = "^ a^^piPj + Wu = 
Sq + VF;7 be the special case that is in involution with [pf +p\+ /\ + U . Then 

S = So ij'^ 

is the corresponding symmetry of H. Since one can always add a constant to a nondegenerate 
potential, it follows that 1/U defines an inverse Stackel transform of H to H. See j37j for many 
examples of this transform. 

9 Mult isepar ability and Stackel equivalence 

From the general theory of variable separation for Hamilton-Jacobi equations jl6| I17j we know 
that second order symmetries Si, S2 define a separable system for the equation 

pI + pI+pI 

H = + V{x, y,z) = E 

X{x,y,z) 

if and only if 

1. The symmetries Ti, Si, S2 form a linearly independent set as quadratic forms. 

2. {5i,52} = 0. 

3. The three quadratic forms have a common eigenbasis of differential forms. 

This last requirement means that, expressed in coordinates x, y, z, at least one of the matrices 
A(^j-^{x) can be diagonalized by conjugacy transforms in a neighborhood of a regular point and 
that [^(2)(^)i -^(i)!^)] = 0- However, for nondegenerate superintegrable potentials in a confor- 
mally flat space we see from Section [S] that 

{5i,52} = [^(2)(a;o),^(i)(a;o)] = 0, and T{xo) = 

so that the intrinsic conditions for the existence of a separable coordinate system are simplified. 

Theorem 7. Let V be a superintegrable nondegenerate potential in a 3D conformally flat space. 
Then V defines a multiseparable system. 

See lU for the details of the proof. 

In [HHI the following result was obtained. 

Theorem 8. Let ui, U2, U3 be an orthogonal separable coordinate system for a 3D conformally 
flat space with metric ds^ Then there is a function f such that 

fds^ = ds^ 

where ds'^ is a constant curvature space metric and ds'^ is orthogonally separable in exactly these 
same coordinates ui, U2, u^. The function f is called a Stackel multiplier with respect to this 
coordinate system. 

Thus the possible separable coordinate systems for a conformally flat space are all obtained, 
modulo a Stackel multiplier, from separable systems on 3D flat space or on the 3-sphere. This 
result provides evidence that, as in the 2D case, all nondegenerate 3D superintegrable systems 
on conformally flat spaces are Stackel equivalent to a superintegrable system on either 3D flat 
space or the 3-sphere, but we have not yet settled this issue. 
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10 Discussion and conclusions 

We have shown that all classical superintegrable systems with nondegenerate potential on real 
or complex 3D conformally flat spaces admit 6 linearly independent second order constants of 
the motion (even though only 5 functionally independent second order constants are assumed) 
and that the spaces of fourth order and sixth order symmetries are spanned by polynomials in 
the second order symmetries. (An interesting issue here is the form of the functional dependence 
relation between the 6 linearly independent symmetries. It appears that the relation is always 
of order 8 in the momenta, but we have as yet no general proof.) This implies that a quadratic 
algebra structure always exists for such systems. We worked out their common structure and 
related it to algebras of 3 x 3 symmetric matrices. We demonstrated that such systems are always 
multiseparable, more precisely they permit separation of variables in at least three orthogonal 
coordinate systems. 

We also studied the Stackel transform, a conformal invertible mapping from a superintegrable 
system on one space to a system on another space. Using prior results from the theory of 
separation of variables on conformally flat spaces, we gave evidence that every nondegenerate 
superintegrable system on such a space is Stackel equivalent to a superintegrable system on 
complex Euclidean space or on the complex 3-sphere, though we haven't yet settled the issue. 
This suggests that to classify all such superintegrable systems we can restrict attention to these 
two constant curvature spaces, and then obtain all other cases via Stackel transforms. We are 
making considerable progress on the classification theory though the problem is complicated. 

All of our 2D and 3D classical results can be extended to quantum systems and the Schrodinger 
equation and we are in the process of writing these up. 

Another interesting set of issues comes from the consideration of 3D superintegrable systems 
with degenerate, but multiparameter, potentials. In some cases such as the extended Kepler- 
Coulomb potential there is no quadratic algebra, whereas in other cases the quadratic algebra 
exists. Understanding the underlying structure of these systems is a major challenge. Finally 
there is the challenge of generalizing the 2D and 3D results to higher dimensions. 
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